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I ABSTRACT 


^Residual  stresses  are  derived  for  a transversely  anisotropic  thick 
hollow  cylinder  which  has  been  chemically  vapor  deposited  at  an  elevated 
ten^erature.  Such  stresses  arise  because  of  the  differential  rates  of 
contraction  in  the  radial  aid  tangential  directions  and  the  anisotropic 
elastic  constants.  Residua  stress  distributions  for  cylinders  with  a 
wall  ratio  Couter  to  inuor  radius)  of  1.30  of  pyrolytic  graphite  and 
pyrolytic  silicon  carbide  (A-SiC)  are  presented  as  a function  of  the 
radius  to  inner  radius.  Thk  effect  of  the  variation  of  the  elastic 
anisotropy  on  the  tangential ^stress  at  the  inner  and  outer  radii  is  pre- 
sented as  a function  of  the  w«l  ratio.  Finally,  the  tangential  and 
axial  stresses  at  the  inner  and  outer  radii  and  the  maximum  radial  stress 
of  chemically  vapor-deposited  oi-SiC  are  presented  as  a function  of  the 
wall  ratio. 
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inner  radius  of  cylinder 
outer  radius  of  cylinder 

elastic  stiffnesses,  i.e.,  cn,  ciz,  C13,  C14,  C33,  C44 
(see  Appendix  C) 

anisotropy  parameter 

radius  coordinate 

elastic  coii|>liances,  i.e.,  sn,  si2»  S13,  S14,  S33,  S44 
(see  Appendix  C) 

thickness  of  thin-walled  cylinder 
displacement  in  the  radial  direction 
wall  ratio  = b/a 
axial  coordinate 

Young's  modulus  in  the  6 direction 
Young's  modulus  in  the  r direction 
Young's  modulus  in  the  z and  6 directions 
ambient  temperature 
deposition  ten^erature 

contractual  strain  in  the  0 direction  due  to  cc’>l-down 

contractual  strain  in  the  r direction  due  to  cool-down 

strain  in  0 direction 

strain  in  the  r direction 

strain  in  z direction 

angular  coordinate 

Poisson's  ratio  in  the  0 direction  due  to  a stress  in 
z direction 

Poisson's  ratio  in  the  z direction  due  to  a stress  in 
0 direction 

Poisson's  ratio  in  the  r direction  due  to  a stress  in 
the  0 direction 

Poisson's  ratio  in  'che  0 direction  due  to  a stress  in 
the  r direction 

Poisson's  ratio  in  the  r direction  due  to  a stress  in 
the  z direction 

Poisson's  ratio  in  the  z direction  due  to  a stress  in 
the  r direction 

stress  in  the  8 direction  (tangential  stress) 
stress  in  the  r direction  (radial  stress) 
stress  in  the  z direction  (axial  stress) 
tangential  stress  ratio 
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I.  INTRODUCTION 


Silicon  carbide  (SiC)  and  silicon  nitride  (Si3N4)  are  candidate  materials 
for  high  temperature  applications  in  ceramic  gas  turbine  for  vehicles^  and  elec- 
tric power  generation.^  The  usual  method  of  manufat. uring  these  materials  into 
intricate  shapes  is  by  hot  pressing  billets  and  diamond  grinding.  Such  a finish- 
ing process  is  difficult  and  expensive  because  of  the  high  hardness  of  these  mate- 
rials. Alternative  methods  of  manufacturing  silicon  carbide  and  silicon  nitride 
have  been  investigated  to  reduce  fabrication  costs.  One  of  these  methods  that 
appears  to  have  potential  promise  is  the  chemical  vapor  deposition  (CVD)  process, 
because  it  readily  allows  the  formation  of  complex  shapes  of  ceramic  materials 
on  preshaped  substrates.  This,  of  course,  minimizes  mechanical  fabrication 
procedures . 

Vapor  deposition  is  the  formation  of  a solid  deposit  occurring  as  condensa- 
tion of  elements  or  compounds  from  the  vapor  state.  Chemically  vapored  deposits 
are  formed  by  chemical  reactions  which  take  place  on,  at,  or  near  the  deposition 
surface,  i.e.,  mandrel  or  substrate.  Deposition  temperatures  of  structural  mate- 
rials of  interest  are  relatively  high;  for  example,  pyrolytic  silicon  carbide 
can  be  formed  at  temperatures  as  high  as  1800  C (3272  F)  and  pyrolytic  silicon 
nitride  up  to  1500  C (2732  F) . The  usual  method  of  producing  pyrolytic  shapes 
of  silicon  carbide  and  silicon  nitride  is  to  introduce  the  appropriate  gas  mixture 
into  a chamber  containing  the  heated  mandrel,  allow  evacuation  of  the  exhaust 
gases,  and  with  proper  control  of  the  complex  physical  and  chemical  steps,  depo- 
sition will  occur.  For  a detailed  description  of  the  CVD  silicon  carbide 
process  see,  for  example,  Weiss ^ and  for  the  CVD  silicon  nitride  process  see 
Niihara.** 

Although  CVD  silicon  carbide  or  silicon  nitride  shapes  are  relatively  easy 
to  produce,  their  tensile  strengths  are  less  than  that  of  respective  hot -pressed 
materials.  This,  in  part,  is  due  to  residual  stresses  inherent  in  the  process 
of  producing  vapor-deposited  materials  which  arise  because  of  several  possible 
mechanisms.^»^  Even  though  physical  processes  such  as  temperature  gradient  and 
structural  growth  or  phase  changes  during  deposition  may  be  mathematically  trac- 
table, the  only  mechanism  considered  here  resulting  in  residual  stress  is  that 
caused  by  the  therma-  anisotropic  coefficients  of  expansion  during  uniform 
cool-down. 

The  complete  detailed  derivation  for  the  thick -walled  cylinder  case  is  given 
in  Appendix  A.  A thick  hollow  cylinder  configuration  was  chosen  as  the  geometry 
to  analyze  because  such  a configuration  is  compatible  with  cylindrical  anisotropy 
and  also  provides  an  ideal  vehicle  by  which  the  residual  strains  and  thus  the 


1 . McLEAN,  A.  F.  Ceramics  In  Small  Vehicular  Gas  Turbines  in  Ceiunin  for  High  Perfoimance  Applications.  Proc  of  the  tad  Army 
Materials  Technology  Conf.,  ed.  by  J.  J.  Burke,  A.  E.  Gorum,  and  R.  N.  Katz,  Brock  mi  Publ.  Co.,  1974,  p.  9-36. 

2.  BRATTON,  R.  J.  Ceramics  in  Gas  Turbines  for  Electrical  Power  Generation  In  Ceramics  for  Hl^  Performance  Apptications.  Proc. 
of  the  2nd  Army  Materials  Technology  Conf.,  ed.  by  i.  i.  Burke,  A.  E.  Gorum,  and  R.  N.  Katz,  Btook  Hill  Publ.  Co.,  1974, 

p.  37-60. 

3.  WEISS,  J.  The  Relationship  of  Structure  ant  Properties  to  Deposition  Conations  and  the  Origin  of  Residual  Stress  in  Chemically 
Vapor  Deposited  Silicon  Carbide.  Ph.D.  Thedt,  Rensselaer  Polytech'^.lc  Inst.,  1974,  p.  49-53. 

4.  NllHARA,  K.,  and  HIRAI,  T.  Chemical  Vapor-Deposited  Silicon  Nttride.  1.  of  Mater.  Sd.,  v.  11, 1976,  p.  593-603. 

5.  POWELL,  C.  F.,  OXLEY,  J.  H..  and  BLOCHER.  I.  M..  Jr.  Vapor  Deposition.  John  Wiley  and  Sons.  Inc.,  1966,  p.  659-662. 


stresses  can  be  experimentally  determined.*  The  resulting  formulae  are  reduced 
so  that  they  are  applicable  to  a thin-walled  cylinder  and  these  equations  are 
given  in  Appendix  B. 

In  order  to  provide  numerical  results,  two  materials  were  considered:  pyro- 
lytic graphite  and  pyrolytic  silicon  carbide.  Pyrolytic  graphite  was  developed 
and  exploited  for  use  in  missiles  during  the  early  1960 's.  Its  properties  are 
well  characterized  and  are  given  in  Appendix  C.  Calculations  lesulting  from  such 
well-defined  properties  will  reflect  the  accuracy  of  the  data  base.  However, 
there  is  a paucity  of  experimental  pyrolytic  silicon  carbide  (a-SiC)  property 
data  as  indicated  in  Appendix  C,  and  idealized  values  of  the  anisotropic  elastic 
constants  were  used  in  determining  the  distribution  and  magnitudes  of  the  resid- 
ual stresses.  Therefore,  the  residual  stress  calculations  are  those  that  would 
result  when  the  material  (o-SiC)  has  high  preferential  orientation  occurring  at 
a deposition  temperature  of  1800  C (3272  F)  and  thus  these  values  are  considered 
as  upper  bounds. 

As  previously  mentioned,  the  detailed  steps  of  mathematical  analysis  is 
given  in  Appendix  A and  will  not  be  repeated  here.  However,  a brief  description 
of  the  physical  problem  and  analysis  is  given  below,  as  well  as  resulting  perti- 
nent equations. 


II.  ANALYSIS 

The  geometry  considered  is  an  infinitely  long,  thick  hollow  cylinder  with 
an  inner  radius  a and  outer  radius  b.  A cylindrical  coordinate  system  is  used 
with  the  three  normal  directions  0,  r,  and  z,  as  shown  in  Figure  1.  It  is  assumed 
that  the  material,  deposited  at  the  elevated  temperature,  is  preferentially 
aligned  in  the  radial  direction.  Thus,  there  is  cylindrical  anisotropy.  It  is 
assumed  that  elastic  parameters  remain  constant  during  deposition  and  subsequent 
uniform  cool-down  and  there  is  no  mechanical  interaction  between  the  substrate 
(mandrel)  and  deposit,  i.e.,  no  adhesion  and/or  no  expansion  coefficient  mismatch. 
Such  problems  can  be  handled  independently,  see  Reference  6 for  example.  Also, 
during  deposition  the  temperature  throughout  the  body  is  constant  and  uniform. 


*As  of  this  writing  no  experimental  lesultt  have  been  obtained. 

6.  RUDDLESEN,  S.  N.  Thermal  Stresses  in  Pyrolytic  Deposits  in  Special  Ceramics  4,  ed.  by  P.  Poppet,  British  Ceramic  Research 
Assoc.,  1968,  p.  153-163. 


The  analysis  given  in  Appendix  A consists  of  formulating  the  stress -strain 
relationships  of  a body  with  cylindrical  anisotropy  and  substituting  the  relation- 
sli  pj;  for  the  two  stresses  099  and  o^r  into  the  governing  equilibrium  equation. 
Tnis  results,  after  some  manipulation,  in  the  formulation  of  an  equidimensional 
Euler  or  Cauchy^  second-order  linear  differential  equation.  This  differential 
equation  is  then  solved  to  yield  radial  displacement  u as  a function  of  the  radius 
r.  Finally,  through  the  use  of  the  stress-strain  relationships  and  the  boundary 
conditions  the  stresses  099,  a^Tt  and  azz  are  derived.  The  equations  for  these 
stresses  are  given  in  the  following: 
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7.  HILDKBRAND,  I'.  B.  Advanced  Calculus  for  Applicattons.  Prentice-Hall.  1962,  p.  13. 


The  tangential,  radial,  and  axial  strains  are  given  by  the  following  formulas.: 
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The  stresses  099  and  Ozz  at  r = a and  r = b were  also  determined  and  are: 
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From  inspection  it  was  determined  that  the  absolute  value  of  maximum  tangen- 
tial and  axial  stress  occurs  at  r = a.  Also,  the  absolute  maximum  radial  stress 
was  obtained  by  maximizing  Equation  2 and  found  to  be  located  at: 
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The  tangential  and  radial  strains  at  the  inner  and  outer  radii  can  also  be 
obtained  and  are: 
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In  Appendix  B the  formulae  describing  residual  tangential  and  axial  stresses 
at  the  inner  and  outer  radii  can  be  reduced  to  those  of  thin-wall  cylinder  theory. 
Their  usage  is  dependent  upon  the  error  one  is  willing  to  tolerate  and  that  error 
is  dependent  on  the  wall  ratio  and  the  material,  in  particular,  the  value  of  m. 
Refer  to  Section  IV  for  a detailed  discussion  of  this.  However,  for  the  sake  of 
completeness  these  equations  are  given  by  the  following: 


(a6e)j._ 
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In  the  next  section  both  the  more  exact  equations  and  the  thin-wall  theory 
equations  are  applied  to  examples  of  cylinders  of  two  different  materials  having 
a wall  ratio  of  1.30  and  a great  difference  in  anisotropic  elastic  constants. 

The  tangential  stress  ratio  oqq/C  at  the  inner  and  outer  radii  is  evaluated  as  a 
function  of  the  wall  ratio  and  for  a wide  range  of  the  anisotropy  parameter  m. 
Finally,  the  residual  tangential  and  axial  stress  at  the  inner  and  outer  radii 
as  well  as  the  maximum  residual  radial  stress  is  determined  for  o-SiC  as  a func- 
tion of  the  wall  ratio. 


III.  RESULTS  AND  DISCUSSION 

The  material  properties  of  two  tra.isversely  anisotropic  materials,  pyrolytic 
graphite  and  pyrolytic  silicon  carbide,  were  utilized  to  demonstrate  the  appli- 
cability of  the  pertinent  equations.  Pyrolytic  graphite  was  chosen  because  its 
properties  are  well  documented. Pyrolytic  silicon  carbide  was  chosen  because 
of  its  attractive  potential  use  as  a high  temperature  structural  material.  In 
Reference  10,  it  is  indicated  that  the  anisotropic  elastic  constants  associated 
with  the  hexagonal  silicon  carbide  system  (o-SiC  6H)  are  not  likely  to  be  differ- 
ent from  other  polytypes.  Not  all  of  the  stiffnesses  for  a-SiC  6H,  which  were 
determined  experimentally  in  Reference  10,  were  obtained,  whereas  ail  of  these 

8 DONADIO,  R.  N.,  and  PAPPAS.  J.  Mechanical  Properties  of  Fyroistu  Graphite  Raytheon  Techniial  Memorandum  1547,  1964 

9.  High  Temperature  Materiab  Inc.,  Data  Sheet,  October  iSfcO. 

10.  ARLT,  G , and  SCHODDER,  G V.  Some  Elasttc  Cons'er  is  of  Silicon  Carbide  1.  Acoust  Soc  of  Am  , v 37,  no  ’.let  ^ry  1%5. 
p.  384-386 


for  the  trigonal  system,  which  were  calculated  from  idealized  consideration,  were 
given  and  compared  closely  to  those  that  were  experimentally  determined  for  a-SiC 
6H.  The  compliances  for  the  trigonal  system  were  computed  via  the  use  of  inter- 
conversion equations  given  in  Reference  11  and  thus  the  anisotropic  elastic  con- 
stants for  a trigonal  rather  than  hexagonal  silicon  carbide  system  vjere  utilized 
in  subsequent  residual  stress  determinations ^ 

The  coefficients  of  expansion  for  hexagonal  silicon  carbide  (a-SiC)  were 
obtained  from  Reference  12  and  are  also  used  in  subsequent  calculations.  The 
details  of  how  the  constants  were  obtained  are  given  in  Appendix  C.  Anisotropic 
elastic  constants  and  other  constants  such  as  Bee,  Brr,  and  m used  in  residual 
stress  calculations  for  pyrolytic  graphite  and  pyrolytic  silicon  carbide  are 
given  in  Table  1,  as  well  as  the  reference  source,  where  applicable. 

It  was  intended  to  examine  pyrolytic  silicon  nitride  as  well.  However,  a 
search  of  the  literature  indicated  that  there  were  no  anisotropic  elastic  con- 
stant data  available. 

The  stress. distributions  for  oee , oj-r,  and  Ozz,  according  to  Equations  1 to 
3,  as  a function  of  r/a  and  a wall  ratio  of  1.3  ^or  both  pyrolytic  graphite  and 
pyrolytic  silicon  carbide  are  shown  in  Figure  2.  Although  pyrolytic  graphite 
cannot  withstand  such  high  self-imposed  stresses  and  would  fracture  at  a much 
smaller  wall  ratio,  the  results  are  presented  to  compare  the  behavior  of  the  two 
materials. 


Table  1.  ANISOTROPIC  ELASTIC  CONSTANTS 


Constant 

Pyrolytic 

Graphite 

Ref. 

Pyrolytic 

Silicon  Carbide 

Ref. 

Eeo 

4.29  X 10®  psi 
(29.58  X 10®  MN/m^) 

8 

60.6  X 10®  psi 
(417.82  X 10®  MN/ir") 

★ 

Err 

1.55  X 10®  psi 
(10.69  X 10®  MN/m^) 

8 

74.1  X 10®  psi 
(510.9  X 10®  K:i/m^) 

* 

'’00 

-0.15 

8 

+0.255 

♦ 

''or 

+0.90 

8 

-rO.079 

★ 

'’ro 

+0.325 

t 

+0.097 

+ 

Soo 

-4600  X 10'®  in. /in. 

T(j  of  2150  C (3902  F) 

9 

-9624  X 10-®  in. /in. 

T(j  of  1800  C (3298  F) 

12 

6rr 

-51,700  X 10-®  in. /in, 
Td  of  2150  C (3902  F) 

9 

-9204  X 10"®  in. /in. 

Td  of  1800  C (3298  F) 

12 

m 

1.415 

t 

0.932 

t 

♦Indirectly  obtained  from  Reference  10 

fCalculated  from  thn  reripr^cal  relationship  vre  - vor  Err/Eo*) 

Tm  = f t u„/[rr(  1 •vyrvptj/l -Van)]  * 


11  lit  VRMON  R I S In /’iiv  )i/u(  n..fi  lo  Applied  AiH^otropu  l-la^neii%  Oxford  U Ftes>.  1961,  p 25 

12  lAyiOK  A in, I lilNfs  k M Ihe  Irwial  Stnuiure  and  Ifiermal  t tpanuon  of  Cubit  and  Hexagonal  Sdu  on  Carbide  Prot  of 
'hi  ( .oit  . o Silk  in  * ifhi'U  I't-Ti’inion  Prvsx  1960  p 14MS4 


It  is  not  to  be  inferred  here  that  the  stress  distributions  shown  in  Figure 
2 are  absolute.  That  is,  mechanisms  such  as  temperature  gradient,  phase  trans- 
formations, for  example,  will  impose  other  residual  stress  patterns.  The  residual 
stress  distributions  presented  here  are  caused  by  only  one  of  the  various  acting 
mechanisms  and  are  a part  of  the  sum  total. 

Notice  that  oqq  and  cjrr  as  a function  of  r/a  shown  in  Figure  2a  are  of  oppo- 
site sign  compared  to  the  corresponding  stresses  in  Figure  2b.  This  is  caused 
by  the  constant  C,  recalling  that 

c - [grr-(l-»-Vflfl-Vflr)3flfl]Eflfl  ^ 

, 2 
i-vee 

The  different  sign  of  the  tangential  and  radial  stresses  for  the  two  materials 
is  due  to  the  relative  magnitudes  of  the  radial  and  tangential  shrinkage  during 
cooling.  That  is,  if  grr  < (l+vge-ver) 300  which  is  typical  of  pyrolytic  graphite, 
and  if  Err  > (1+V00-V0r)300  which  was  determined  for  pyrolytic  silicon  carbide, 
then  the  tangential  and  radial  stress  distributions  will  be  of  opposite  sign. 

The  axial  stress  distributions  for  the  two  materials  are  of  the  same  sign  because 
V00  of  the  two  materials  is  of  opposite  sign,  thus  compensating  for  opposing  sign 
of  the  constant  C associated  with  each  material. 

The  absolute  values  of  the  stress  magnitudes  for  pyrolytic  graphite  are 
greater  than  those  of  pyrolytic  silicon  carbide  for  the  same  wall  ratio.  This 
is  mainly  attributed  to  the  constant  C,  which  for  pyrolytic  graphite  is  very  much 
greater  than  that  for  pyrolytic  silicon  carbide.  Again,  the  magnitude  of  C is 
strongly  dependent  upon  the  contraction  rates  during  cooling,  i.e.,  3rr  $00, 
as  well  as  E©©.  Since  6rr  for  pyrolytic  graphite  is  much  greater  than  that  of 
pyrolytic  silicon  carbide,  then  pyrolytic  graphite  will  have  greater  stress 
magnitudes  even  though  £©9  is  an  order  of  magnitude  less  than  that  of  pyrolytic 
silicon  carbide. 

It  is  interesting  to  note  than  even  if  Brr  = 300,  residual  stresses  would 
still  persist  because  of  the  elastic  anisotropy. 

The  influence  of  the  elastic  anisotropy  parameter  m on  the  tangential  stress 
at  the  inner  and  outer  radii  of  a thick  cylinder.  Equations  7 and  8,  was  evaluated 
as  a function  of  the  wall  ratio  W and  is  shown  by  solid  curves  in  Figure  3.  Prac- 
tical values  of  m ranging  from  approximately  0.8  to  1.6  for  hexagonal  and  trigonal 
crystal  structures  were  obtained  from  References  11  and  13.  Therefore,  m was 
varied  from  0.50  to  2.0  as  shown  in  Figure  3.  Note  tliat  beyond  a wail  ratio  of 
approximately  1.30  the  absolute  magnitude  of  the  tangential  stress  ratio  increases 
as  m decreases.  At  wall  ratios  of  approximately  less  than  1.30  the  magnitude  of 
m appears  to  have  little  effect  on  the  tangential  stress  at  the  inner  and  outer 
radii . 


13  HtARMON,  R F S The  Elastic  Constants  of  Anisotropic  Materials  ■ If  Advances  in  Physics  - A Quarterly  Supplement  of  the 
Phil.  Mag  , V 5,  no  19,  July  1956,  p.  323-381. 
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Figure  2.  Residual  stress  distribution  as  a function  of  r/a  for 
anisotropic  thick  hollow  cylinders  of  (a)  pyrolytic  graphite 
and  (b)  pyrolytic  silicon  carbide. 
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Figure  3.  Tangentiai  stress  ratio  at  inner  and  outer  radius  as  a 
function  of  the  waii  ratio  W,  and  the  anisotropy  parameter  m. 

In  Reference  14  the  residual  stresses  in  transversely  anisotropic  hexaferrite 
material  during  cooling  were  determined  ignoring  Poisson's  ratios  and  axial 
deformation.  This,  of  course,  simplifies  the  analysis  but  dependent  upon  the 
material  constants,  the  error  could  be  extreme.  For  example,  for  the  special 
case  V00  = V0r  = 0,  m = (E00/Err)^>  C = (err-380)E00,  and  applying  this  approach, 
the  stresses  would  differ  from  values  calculated  from  more  exact  theory  by  -50% 
on  a cylinder  of  pyrolytic  graphite  and  -86%  on  one  of  pyrolytic  silicon  carbide. 

Figure  4 presents  the  tangential  and  axial  stresses  at  the  inner  and  outer 
radii  as  well  as  the  maximum  radial  stress  as  a function  of  the  wall  ratio  for 
pyrolytic  silicon  carbide  deposited  at  1800  C (3272  F) . The  maximum  radial 


14.  KOOLS,  F.  Complications  in  Firing  Oriented  Hexaferrites  Due  to  Anisotropic  Effects  Cracking  of  Radially  Oriented  FHnges  During 
Firing.  Science  of  Ceramics,  published  by  the  Societe  Prancaise  De  Ceramique,  v.  7,  1973,  p.  27-45. 


Figure  4.  Residual  tangential  and  axial  stresses  at  the  inner 
and  outer  radii  and  maximum  radial  stress  for  a thick 
hollow  cylinder  of  pyrolytic  silicon  carbide  as  a function 
of  the  wall  ratio  W. 


stress  was  determined  by  first  obtaining  the  location  of  this  stress,  from  Equa- 
tion li  and  using  that  value  in  Equation  2.  The  reader  is  cautioned  that  the 
stresses  shown  represent  upper  bounds  in  that  it  is  assumed  that  the  material  is 
fully  anisotropic.  It  appears  from  Figure  4 that  the  largest  tensile  stres.’  will 
occur  in  the  tangential  direction  at  the  outer  radius  of  the  cylinder.  If  the 
mechanism  of  differential -direction  cooling  due  to  anisotropy  of  C\T)  silicon  car- 
bide is  the  major  cause  of  residual  stress,  then  a longitud.-.nal  fracture  will 
occur  and  originate  at  the  outer  radius.  Thus,  for  a given  material  and  deposi- 
tion conditions  there  should  result  a critical  wall  ratio  above  which  fracture 
will  occur,  see  Reference  14. 

In  Section  II  the  tangential  and  axial  stress  equations  were  siraplified 
assuming  the  cylinder  wall  was  thin,  resulting  in  Equations  16  and  17.  The  stress 


ratios  (O09/C)j.=a  (o00/C)y-|,  based  on  Equations  16  and  17  for  m values  of 

0.50  and  2.00  are  shown  plotted  as  dashed  lines  in  Figure  3 to  facilitate  compar- 
ison to  the  more  exact  theory  given  by  the  solid  curves.  It  is  readily  seen  that 
the  differences  between  the  more  accurate  thick-wall  cylinder  theory  (solid  curves) 
and  the  thin-wall  cylinder  theory  (dashed  lines)  results  can  be  quite  extreme  and 
depend  upon  the  wall  ratio  W and  the  anisotropy  parameter  m. 


IV.  SUMMARY 

1.  The  signs  of  the  stresses  are  dependent  upon  the  magnitudes  of  the 
directional  contractions  during  cooling,  and  300,  and  Poisson's  ratios,  V00 
and  V0r. 


2.  Even  if  the  contraction  during  cooling  in  the  radial  (3rr)  an*!  tangential 
(300)  directions  were  equal,  residual  stresses  would  persist  because  of  elastic 
anisotropy. 

3.  The  absolute  maximum  stress  which  results  from  uniform  cool-down  is  the 
tangential  stress  at  the  inner  radius  of  the  cylinder  for  all  materials  of  trans- 
verse (orthotropic)  anisotropy. 

4.  The  maximum  tensile  stress  for  a fully  anisotropic  cylinder  of  pyrolytic 
silicon  carbide  material  (o-SiC)  occurs  in  the  tangential  direction  at  the  outer 
radius.  Fracture  should  originate  at  this  location  and  extend  in  the  longitudinal 
and  radial  directions.  This  mechanism  should  result  in  a critical  wall  ratio 
above  which  fracture  will  occur. 

5.  As  the  elastic  anisotropic  parameter  m decreases  the  tangential  stress 
increases  when  W > 1.30.  If  W is  less  than  1.30,  m will  have  little  effect  on 
the  magnitude  of  this  stress. 

6.  If  axial  deformation  and  Poisson's  ratios,  V00  and  V0x,  are  ignored, 
errors  in  determining  the  residual  stresses  are  dependent  upon  the  properties  of 
the  material.  The  stresses  would  differ  from  those  values  calculated  from  more 
exact  theory  by  approximately  -50%  for  pyrolytic  graphite  and  -86%  for  pyrolytic 
silicon  carbide. 

7.  The  differences  between  the  thin-wall  cylinder  theory  and  the  more 
accurate  thick-wall  cylinder  theory  results  can  be  quite  extreme  and  depend  upon 
the  wall  ratio  W and  the  anisotropic  parameter  m. 


APPENDIX  A.  ANALYSIS 


The  geometry  considered,  as  previously  mentioned,  is  a long,  thick  hollow 
cylinder,  as  shown  in  Figure  1,  with  the  three  direction  coordinates  6,  r,  and  z. 
The  hexagonal  crystal  structure  is  assumed  to  deposit  at  an  elevated  ten^erature 
in  a preferred  orientation  in  the  radial  direction  such  that  there  is  isotropic 
symmetry  about  the  Z-axis  and  anisotropy  exists  in  the  radial  direction.  It  is 
also  assumed  that  the  elastic  constants  remain  constant  during  deposition  and 
subsequent  cool-down  is  gradual  such  that  the  temperature  throughout  the  cylinder 
is  uniform.  It  is  further  assumed  that  there  is  no  mechanical  interaction  between 
the  substrate  (mandrel)  and  deposit,  i.e.,  no  adhesion  or  no  expansion  coefficient 
mismatch.  Also  during  deposition  there  is  no  temperature  gradient  occurring  within 
the  body  (this  problem  can  be  solved  separately) . 


According  to  Reference  15  the  stress-strain  relationships  of  such  a body  with 
cylindrical  anisotropy  is: 


€00  = (l/E00)(o00-vzecTzz)-(vr9/Err)cfrr 
Err  = (l/Err)®rr-(\’er/E00)cr00-(vzr/E00)ozz  ^ 
Ezz  * (VEzz)ozz”(''0z/^00)'^00"('’rz/Err)orr 


(A-a) 


Notice,  because  of  the  conditions  of  symmetry,  there  are  neither  shear  stresses 
nor  shear  strains.  Also,  since  Ezz  * Eee,  vze  * vez  = V00,  = vrz»  '’zr  “ '’0r» 

and  from  the  theorem  of  reciprocity  we  have  - '’0r/E00»  thus: 


E00 

Err 


(I/E00)  (a00-vo6Ozz-ver'^rr) 
(1-Err)  [orr'^’rO  (®00'’'®zz)  ] 1 


Ezz 


(I-E00)  (a2,z-V00a00-V0rOrr) 


(A-b) 


During  uniform  cool-down  the  contraction  of  the  cylinder  in  the  r direction 
is  different  than  that  either  in  the  0 or  z direction  because  of  the  differing 
thermal  coefficients  of  expansion  of  the  material.  Thus,  additional  strain  com- 
ponents must  be  added  to  Equation  A-b  as  given  in  the  following: 


'=00  = 
Err  = 


(I/E00)  (a00-V00azz-ver«Jrr)+P00 

(1/Err)  [orr"'’r0(®00'‘'®zz)]'*‘Brr 


Ezz 


(I/E00)  (ozz-V00O00-V0i.Orr)+E00 


(A-c) 


where  $00  and  Srr  is  the  total  contraction  in  the  tangential  and  radial  direction 
due  to  isothermal  cool-down. 


Ne  shall  first  let  Ejz  = 0 and  later  modify  the  solution  to  con^jly  with  an 
infinitely  long  cylinder  with  a stress-free  end  condition,  thus  Equation  A-c 
becomes : 


15.  LEDHNnSXII,  S,  G.  Theory  of  EUuHcity  of  an  Anisotropic  Elastic  Body  In  Holden-Day  Scries  in  Mathematicel  Physics,  ed.  by 
J.  L.  Biandstatter,  1963,  p.  25. 


<^0e-'>0r<^rr-V0e(ozz)e^^,O  * ^98 
Vr0<T00-arr+Vr0(ozz)c__«o  * 

xz 

V00a00+verar-(ozz)e22«O  * ^ 

where 


(A-d) 


^00  * (e00-P00)E00 

^rr  * “(Err-Prr)^!^:* 
^ » 300E00 


Equation  A-d  can  be  solved  such  that  the  normal  stresses  oqq,  ^d  Ozz 
can  be  expressed  in  explicit  form;  when  this  is  accomplished  we  have: 


000  = fl/D)[(l-V0rVr0)X00-ver(l+v00)Xrr-(ve0+v0rvr6)Xl 
Orr  » (l+vea)/d[vr8A00-(l-V00)An-vr0X] 

Uzz)e22=0  “ (l/D)[(''00+vr0V0r)X00-(l+vO0)v0rXrr-{l-Vr0'>0r)X]^ 


(A-e) 


where 


D * (l+V00)(l-V00-2v0rVi.0) 

The  equilibrium  equation  for  deformation  symmetrical  about  the  axis  of  a thick- 
walled  cylinder^®  is: 


dai.j./dr+(ai.i.-O00)/r  = 0 (A-f) 

and  also  by  definition: 

£00  = u/r,  and 
Ej-y  = du/dr 

It  is  assumed  that  the  constants  B00,  Brr,  E00,  and  Err  are  not  functions  of  the 
radius.  Therefore,  by  differentiating  Orr  with  respect  to  r from  the  second  of 
Equations  A-e  and  substituting  that  result  as  well  as  the  relationships  for  000 
and  Oyy  from  Equation  A-d  into  Equation  A-f  gives: 

r^du/dr^+rdu/dr-ra^u  = a^r  (A-1) 

where 

= (E00/Err)  (l-V0r'’r0)/(l-ve0),  and 

ai  = (1/1-V00) [(l-V90-V0r)8rr-(l-2vr0) (E00/Err)306] 

Equation  A-1  is  an  equidimensional  linear  differential  equation,  which  is  variously 
called  Euler’s  equation  or  Cauchy’s  equation,’^  whose  homogeneous  solution  is; 


16.  TIMOSHENO,  S..  and  0<X)D1ER,  J.  N.  TTiiory  of  Elatticity,  "Ixti  ed.  McGraw-Hill  Book  Co.,  Inc.,  1951,  p.  236. 


U||  ■ Cir®+C2/i*, 
and  the  particular  solution  is 


Up  = air/(l-m2),  m / 1.0. 


Also  note  that  in  t!;,e  above  solution  it  is  assuiced  that  V0jVr0  < 1.0.  The  complete 
solution  is: 

u = Cir®+C2/r*’^+air/(l-m^),  m ^ 1.0.  CA-2) 

Cl  and  C2  are  obtained  from  the  boundary  conditions,  i.e.,  ajT  = 0 when  r = a and 
b.  The  constants  Ci  and  C2  are  determined  through  the  use  of  the  middle  equation 
OjY  of  Equation  A-e,  Equation  A-2,  the  appropriate  strain  definitions,  and  the 
boundary  conditions.  When  this  is  acconqplished  one  determines  the  following 
equations  for  Ci  and  C2: 

C[(l-V00)EQ0/Err-2ver]  (W”*-^-!)^"” 

" E00(l-m2)[vgr+(l-V00)m](l-W=^®) 


and 


C[(l-V0e)Ee6/Err-2v9r3 (W”-^-l)b"^^ 

E00(l-m2)[v0r+(l-V00)m](l-w2m) 


where 

C = [Orr-(l'*’'’00-'>er)Ee9]E00/(l-ve6) 


Substitution  of  Ci  and  C2  into  Equation  A-2  and  dividing  by  the  variable  r, 

gives 

000  = u/r  = C[(l-V00)E00/Err-2V0|]  C 

E00(lV)(l-W^"‘)  \v0r+(l-V00)m  * (A-3) 

(W“-^-l)  W®+1  1 

+ SIT 

[v0r-(l-V00)m]  (r/a)®'*'^  j 

and  by  differentiating  Equation  A-2  with  respect  to  r,  and  substituting  for  Ci 
and  C2  from  the  above,  the  radial  strain  err  can  be  obtained,  which  is: 

err  = du/dr  = mC[(l-V00)E00/Err-2v0|]  f (W®*1-1) 

E00(l-m2)(l-W^nr)  | [V0r+(1-V00)m]  (A-4) 

(W®-^-l)  1 , 

Wai/(i-in2) 

[V0r-(1-V00)ra]  (r/a)®^^J 


Through  the  use  of  Equation  A-e  and  the  definitions  of  strains,  the  tangential, 
radial,  and  axial  stresses  (with  = 0)  can  be  obtained  and  are: 

000  = C/(l-m2)  {l+m/(l-W2”»)[(W®^l-l)(r/a)“"^-(W“”^-l)W®^VCr/a)®*^]}  (A-5) 

Oj^  = C/Cl-m2)  {l+l/(l-w2m)[C#+l-l)(r/a)®-l+(W®-l-l)W®+V(r/a)®*l]}  (A-6) 

(az2)e  =0  ^ {CveeM90+vr0Mrr)Cir®-V(l-m^)  + 

z z 

[('’eeN00+Vr0Ni.r)/(l-m^)]C2/r®'^^+(v00P00+Vr0Pi.r)O00}-!300E00 

where 

M00.=  (l-V0j.Vr0)/(l+VQ0)+mvj0,  N00  = (l-v0rvre)/(l+v00)-myr0 

Mrr  = ver+Cl-^ee)®.  N-rr  = ver-Ci-^ee)® 

Pee  = [(l-V0rVre)/(l+vee)+Vrelai/(l-in^)-(300+Vr0erj.),  and 
Prr  = [C/E00(l-m2)] [(l-v00)E00/Err-2v0j]. 

The  axial  stress  Ozz  can  now  be  determined  such  that  the  stress  at  the  ends  of  the 
cylinder  vanish.  This  is  accomplished  by  obtaining  a resulting  force  C3  in  the 
following  manner: 

C3ir(b2-a2)  = -^(azz)e  irrdr 
J *"zz=0 
a 


Substitution  of  Cazz)ezz=0  Equation  A-7a  into  the  above,  integrating  and  solv- 
ing for  C3  results  in: 


■-2(l+ve0)1J'p00M00+Vr0Mrr1  (W™*1-1)  Cj 

D J\L  (1-ro^)  Cl+m)  J (W^-l)  a^'”* 

veeNee+vreNrrI  (w^~”*-l)  ^2 

^ (l-m2)(l-m)  J (W2-1)  b^'^® 


+ 


Eee 

2 


^v00P00+VrePrr- 


gee  D 

1+V00 


(A- 7b) 


Now  by  adding  C3,  given  by  Equation  A-7b,  to  (Ozz)e2z=0>  given  by  Equation  A-7a 
and  noting  that  M00/M1.1.  = mErr/EGe,  and  N00/N1.1.  = -inEi.i./E0e, 


we  obtain  Ozz,  which  is: 


|Cver*mve9)(W«*>-l)[(r/a)->  - ^ 


2 (W 


W^-l)  J 


(A- 7) 


(v9r-mv06)(W®-l-l)W'”+l 


r—  - 


2 

(1-m)  (W2-1)  _ 


Now  the  axial  strain  Ciz  can  be  obtained  by  substitution  of  Equations  A-5,  A-6, 
and  A-7  into  the  last  of  Equation  A-c.  This  results  in: 


(V6e+V0j.)  2C  f/ver+mv00 

Ee0(l-m2)  " E00(l-m2)(l-w2in)  (W^-l)  I 1+m 


(A- 8) 


Equations  A-3  and  A-4  give  the  tangential  and  radial  strains  as  a function  of 
the  ladius  ratio  parameter  r/a,  the  wall  ratio,  and  the  various  elastic  and 
physical  constants.  The  axial  strain  given  by  Equation  A- 8 is  constant  and 
dependent  only  upon  the  wall  ratio  and  the  elastic  and  physical  parameters. 

Equations  A-5,  A-6,  and  A-7  gives  the  tangential,  radial,  and  axial  stress 
000,  Oj-x,  and  a^z  as  a function  of  the  radius  ratio  parameter  r/a,  the  wall  ratio 
W,  and  the  various  elastic  and  physical  constants. 

The  stresses  and  strains  at  the  inner  and  outer  radii  can  readily  be  obtained 
by  allowing  r = a and  r = b in  the  appropriate  equations.  The  tangential  and 
axial  stresses  at  the  inner  and  outer  radii  are  obtained  from  Equations  A-5  and 
A-7  and  are: 


^(v0x-mv00)(W">-l-l)W"‘-l  jl’ 


(A-9) 


(A- 10) 


(l-m2)(l-w2ra) 


(A-11) 


and 


(®zz) 


(l-m2)  (1-W2®) 


(A-12) 


The  tangential  and  radial  strains  at  the  inner  and  outer  radii  obtained  from 
Equations  A- 3 and  A-4  are: 


(eee) 


and 


C[(l-V99)E99/Err-2v9r]  f 

(W®+^-l) 

(W°-^-l)W®'^^^ 

ai 

r=a 

E09(l-W2®)(l.m2)  \ 

V9r+Cl-V00)m 

V0r-(1-V00)m  J 

1-m^, 

C[(l-V99)E99/Err-2v9|l  ( 

■(W®+1-1)W®"1 

(K®-1_1) 

ai 

r=b 

E9e(l-W2®)(l-m2)  \ 

V0j.+  (1-V00)m 

V00“(l-V00)mj 

, 

1-m^, 

_ mC[(l-V99)E0e/Err-2veJ] 

r a«®^i-i) 

(W®-1-1)W®'"1' 

r=a 

E0e(l-W2“)(i-m'^) 

\v0r+(l-V00)m 

V0r-Cl-V00)m 

n»C[(l-V99)E00/Err-2v9r] 

|'(y^m+l_i)(ym-l 

(W®“1-1) 

1, 

r==b 

' E9e(l.W5^®)(l-m^) 

\'>9i.+  (l-V90)m 

V9r-(1-V00)m 

j l-m2 

(A- 13) 

(A-14) 

tA-15) 


(A- 16) 


Since  e^z  as  given  by  Equation  A-7  is  constant,  there  is  no  difference  in  axial 
strain  at  the  inner  and  outer  radii. 

Examination  of  the  stress  equations  for  099  and  Ozz  reveals  that  the  absolute 
maximum  stress  occurs  at  the  inner  radius,  i.e.,  when  r =«  a. 

The  location  at  which  the  radial  stress  is  a maximum  can  be  found  by  maxi- 
mizing Equation  A-6.  If  this  is  accomplished  we  obtain: 


r/a 


r/m-t-l\  (W»”-^-l) 

tU-i/ 


jyffl+ 1 


hm 


(A-17) 


APPENDIX  B.  THIN  WALL  STRESSES 

The  residual  stresses  009  and  o^z  at  the  inner  and  outer  radii  of  a thin- 
walled  cylinder  can  be  obtained  from  Equations  A-9  and  A-10,  and  A-11  and  A-12. 
This  is  accomplished  by  letting  W = 1+t/a  where  t is  the  cylinder  wall  thickness 
and  using  series  expansions  for  each  of  the  terms  that  appear  in  Equations  A-9  to 
A-12,  For  example,  the  reduction  of  Equation  A-9  to  thin-wall  cylinder  theory  is 
accomplished  in  the  following  manner: 

Note: 

(009)^,  = (C/i-m2)[l+(m/l-W2'“)(2W®+l-l-W2®)],  and 

by  using  a series  expansion  for  and  i.e.: 

= (l+t/a)2n>  = l+2mt/a+m(2m-l) (t/a)2+ (t/a)2  < 1,  and 

= (l+t/a)"**!  = l+(m+l)t/a+m(m+l)(t/a)2/2  + (t/a)2  < l, 

and  substituting  the  appropriate  series  expansion  (only  up  to  the  second-order 
terms)  into  the  above  equation 'as  shown  below  gives: 

^ ^ ®[2|(l+(m+l)t/a+ra(m+l)/2(t/a)2|-l-|[l+2mt/a+m(2m-l 

(’^90)r-a  \T-^  ^ 2m (t/ a)  +m (2m- 1)  (t/ a) ^ 

The  above  reduces  to: 

(O90)j.^g  = -C(t/a)/[2+(2m-l)t/a],  (B-1) 

Utilizing  the  same  approach  the  remaining  desired  equations  are: 

(<^00)j.=j,  = C(t/a)/[2+(2m-l)t/a],  (B-2) 

(ozz)r^a  -V00C(t/a)/2(l+mt/a),  (B-3) 

and  finally: 

(ozz)j.^b  = v0eC(t/a)/2(l+mt/a)  (B-4) 


APPENDIX  C.  ELASTIC  AND  PHYSICAL  CONSTANTS 


The  elastic  constants  for  two  transversely  anisotropic  materials,  pyrolytic 
graphite  and  pyrolytic  silicon  carbiae,  were  used  in  the  body  of  the  report  to 
obtain  numerical  results.  Pyrolytic  graphite  was  chosen  because  its  property  data 
are  well  documented,  see,  for  example.  References  5,  8,  and  9.  Pyrolytic  silicon 
carbide  was  also  considered  because  it  is  a material  of  potential  promise.  Unfor- 
tunately, a search  of  the  literature  revealed  that  there  were  no  anisotropic  prop- 
erty data  available  for  the  other  material  of  interest,  pyrolytic  silicon  nitride. 
Even  though  the  pyrolytic  silicon  carbide  data  from  the  literature  were  incomplete 
the  anisotropic  elastic  constants  were  estimated  indirectly  from  Reference  10. 

The  data  for  pyrolytic  graphite  data  are  first  presented  and  then  that  data 
for  pyrolytic  silicon  carbide  follow. 

The  anisotropic  elastic  constants  utilized  for  calculations  were  those  ob- 
tained frora  Reference  8 as: 

Eee  = 4.29  X 10^  psi  (29,58  x 10^  MN/m^),  vee  = -0.15 

Ej-r  = 1.55  X 10°  psi  (10.69  x 10^  ^'N/m^),  « +0.90 

Although  the  published  value  for  Vj.0  is  0.35,  for  consistency  the  reciprocating 
relationship  was  used,  resulting  in  v^e  = 0.325.  This  value  was  used  throughout 
the  calculations.  Recalling  that: 

ra  = [(Eee/Err)  (l-''0r^r0)/(l-V00) 

and  substituting  the  above  constants  intc  this  I'elationship  gives; 

m = 1.415 

The  total  contraction  during  cooling  in  the  tangential  and  radial  directions 
was  determined  from  Reference  9 as 

000  = 01+02 (Tg-Td) 


where 


ai  = -270  x 10“®  in. /in. 


and 


02  = 1.13  X 10"^  in./(in.-®F) , 


thus 


009  = -27J  X 10"^  in. /in. + [1.13  x 10'^  in./ (i:i. -°F) ] (70  -3902  F) , or 
000  = -4600  X 10"^  in. /in. 


also 


3rr  " 


where 


03  = 13.50  X 10"®  in./(in.-°F), 


and 


Brr  = [13.50  x 10"®  in./ (in. -“F)] (70  F-3902  F) , or 

3j.j.  = -51,700  X 10"®  in. /in. 

The  manner  by  which  the  elastic  and  physical  constants  for  silicon  carbide 
(CVD)  were  calculated  are  now  presented  in  the  following. 

The  authors  of  Reference  10  determined  by  resonance  and  the  double-pulse 
echo  method  for  the  hexagonal  crystal  structure  (o-SiC  6H  (33))  the  magnitude  of 
the  stiffnesses  cu,  ci2i  C44>  C33>  and  cgs  and  the  compliances  su,  si2»  S44  and 
See.  Unfortunately,  no  values  were  reported  for  C13,  S13,  and  S33.  However, 
according  to  Reference  11  there  is  an  interconversion  between  the  stiffnesses  and 
the  compliances  of  a hexagonal  system.  Through  the  use  of  three  simultaneous 
equations,  the  three  unknowns  C13,  sjs,  and  S33  were  calculated.  Using  nominal 
experimental  »:alues  of  the  known  parameters,  S13  was  found  to  be  an  imaginary 
number,  which  can  not  physically  occur.  Closer  examination  of  the  errors  associ- 
ated with  the  experimentally  determined  constants  published  in  Reference  10  re- 
vealed that  extremely  large  errors  would  result  in  the  calculation  of  C13,  S13, 
and  S33  when  using  the  interconversion  equations.  However,  the  elastic  constants 
of  other  polytypes  of  SiC  are  not  likely  t^,  be  very  different  and  have  included 
calculated  values  of  the  transformed  elastic  constants  of  a cubic  polytype  6-SiC.^® 
Although  this  system  is  trigonal  rather  than  hexagonal,  mean  values  of  the  trans- 
formed stiffnesses  compare  quite  well  to  those  experimentally  determined  stiff- 
nesses for  the  hexagonal  structure.  These  transformed  values  of  the  stiffnesses 
for  the  trigonal  system  were  used  here  in  conjunction  with  the  interconversion 
equations  for  a trigonal  system^*  to  determine  the  compliances.  These  calculations 
and  the  results  are  suiranarized  below. 


The 

interconversion  equations 

are: 

(a) 

C11+C12  = S33/X1 

(e) 

C’l 

-C12  ■=  S44/X2 

(b) 

Cl3  = -'^13/^1 

(f) 

= -S14/X2 

(c) 

C33  = (Sll+Si2)/Xi 

(g) 

C44 

= (Sll-Si2)/X2 

(d) 

Xi  = s 13 (si i + Si 2) -2si 3 

(h) 

X2 

= S44 (Si 1 -Si  2 )-2s  14 

Reworking  these  eouations  so  that  they 


Si  . - 

-Ic 

1 3/ U‘  1 i + c 1.  ) ]s  , j 

3 ■ 

!/ 

[c  1 1 v fc  ; , n 

i -- 

(4  i-.  t t.  ; 1 . 

•re  amenable  for  use,  we  obtain: 

7 

Si-  = -Si  *’3^) 

'!l  = - s;  V n/c  n- 3 3) 1 4/ U/<-'44-S44  ) 

s i I.  - k'  i 1 ! -s 


The  transformed  stiffnesses,^®  given  in  lO^^dyn/cm^,  are: 


cii  = 4.793 
Ci2  - 0,981 
Ci3  = 0.558 


Ci4  = 0.598 
C33  = 5.214 
Ci^ “ 1 # 483 


By  utilizing  the  above  equations  the  corresponding  compliances  given  in  units  of 
10" dynes,  are: 


Sji  = 2.392  S24  = 1.211 

$12  “ -0,610  S33  = 1.958 

Si3  = -0.189  S44  = 7.72 

The  remaining  elastic  constants  were  determined  from  the  above  and  are  given 
below: 

E09  = 1/sii  = 60.6  X 10®  psi  (417.82  x lO^  MN/m^) 

Err  = I/S33  = 74.1  X 10®  psi  (510.90  x 10®  MN/m^) 

^00  = V12  = -S12/S11  = +0.255 
V0r  " ^13  = “S13/S11  = +0.079 

and 

^'ro  “ V0r  " Err/E00  = 0.097 

Finally,  m was  deteimined  from  the  above  values  to  be  0.932. 

Reference  12  gives  the  thermal  expansion  curves  in  the  basal  direction  "a" 
and  in  the  "c"  direction  as  a function  of  temperature  for  the  hexagonal  silicon 
carbide  system  (mod  II).  These  curves  were  used  to  determine  the  total  shrinkage 
in  the  0 and  r directions  of  SiC  from  1800  C to  room  temperature  and  were  found 
to  be 

3q0  = -9624  X 10"®  in. /in., 

8j.j.  = -9204  X 10"®  X 10"®  in. /in. 
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National  Aeronautics  and  Space  Administration,  Langley 
Research  Center,  Hampton,  Virginia  2366S 
1 ATTN;  Hr.  H.  F.  Hardrath,  Mall  Stop  188M 
1 Mr.  R.  Foye,  Mall  Stop  I88A 

National  Aeronautics  and  Space  Administration,  Lewis 
Research  Center,  21000  Brookpart  Road,  Cleveland,  Ohio  44135 
I ATTN;  Mr.  S.  S.  Hanson 

I Or.  J.  E.  Srawley,  Mall  Stop  105-1 

1 Hr.  W.  F.  Brown,  Jr. 

1 Mr.  H.  H.  HIrschberg,  Head,  Fatigue  Research  Section, 

Hall  Stop  49-1 

1 Mr.  G.  Mervin  Ault,  Assista  'If'  MAS  Division 

1 Or.  H.  B.  Probst,  MS  49-1 

Office,  Olrector  of  Research  and  De,  .Ofmmnt,  Department 
of  the  Air  Force,  The  Pentagon,  Washington,  D.  C,  20330 
1 ATTN-  AFORO-OR,  LTC  horas  C.  Hamlin 
1 Major  Donald  Sponberg 

National  Bureau  of  Standards,  U.  S.  Department  of  Conaerce, 
Washington,  0.  C.  20234 
I ATTN  Hr.  J.  A.  Bennett 

I Mechanical  Properties  Data  Center,  Belfour  Stolen  Inc., 

13917  W.  Bay  Shore  Drive,  Traverse  City,  Michigan  49634 


No.  of 

Copies To 

Deftnst  Materials  Service,  General  Services  Administration, 
Washington,  0.  C.  20405 

1 ATTN:  Mr.  Clarence  A.  Fredell,  Olrector,  Technical  RAD  Staff 

1 Mr.  W.  F.  Anderson,  Atomics  International,  Canoga  Park, 
California  91303 

Midwest  Research  Institute,  425  Coker  Boulevard,  Kansas  City, 
Missouri  64110 
1 ATTN:  Mr.  G.  Gross 

1 Mr.  A.  Hurileh,  Convair  Division,  General  Dynamics  Corp., 

Hall  Zone  630-01,  P.  0.  Box  80847,  San  Diego,  California  92138 

1 Mr.  J.  G.  Kaufman,  Alcoa  Research  Laboratories,  New  Kensington, 
Pennsylvania  15068 

1 Mr.  P.  N.  Randall,  TRW  Systems  Group  - 0-1/2210,  One  Space  Park, 
Redondo  Beach,  California  90278 
1 Mr.  W.  A.  Van  der  Sluys,  Research  Center,  Babcock  and  Wilcox, 
Alliance,  Ohio  44601 

1 Mr.  B.  M.  Wundt,  2346  Shlrl  Lane,  Schenectady,  New  York  12309 
1 Dr.  Raymond  Bisplinghoff,  Deputy  Administrator,  National 
Science  Foundation,  1800  G.  Street,  NW,, 

Washington,  0.  C.  20550 

1 Dr.  A.  Lovelace,  Deputy  Assistant  Secretary  (RAD),  Office  of 
Assistant  Secretary  of  the  Air  Force  (Research  and  Development), 
Room  4E973,  Pentagon,  Washington,  0.  C.  20330 

1 Or.  A.  R.  C.  Westwood,  RIAS  Division  of  the  Martin  Company, 
Baltimore,  Maryland 

1 Dr,  S.  Wlederhom,  National  Bureau  of  Standards, 

Gaithersburgh,  Maryland  20760 

1 Dr.  P.  Jorgensen,  Stanford  Research  Institute,  Menlo  Park, 
California  34025 

I Mr,  M.  Allen  Hagid,  Materials  Marketing  Engineer, 

Florida  RAO  Center,  Pratt  and  Whitney  Aircraft,  P.  0.  Bax  2691, 
West  Palm  Beach,  Florida  33402 

1 Dr,  E.  A.  Stelgerwald,  TRW  Metals  Division,  P.  0.  Bax  250, 
Minerva,  Ohio  44657 

1 Hr.  H.  J.  Manjoln,  Westinghouse  RAD  Center,  1310  Beulah  Road, 
Pittsburgh,  Pennsylvania  15235 

I Hr,  E.  T.  Wessel,  Westinghouse  RAD  Center,  1310  Beulah  Road, 
Pittsburgh,  Pennsylvania  15235 

1 Or.  Alan  S.  Tntelman,  Failure  Analysis  Associates, 

11777  Mississippi  Avenue,  Suite  4,  Los  Angeles, 

California  90025 

Battelle  Colimibus  Laboratories,  505  King  Avenue, 

Colianbus,  Ohio  43201 
1 ATTN:  Mr.  J.  Campbell 

1 Or.  G.  T.  Hahn 

1 R.  G.  Hoagland,  Metal  Science  Group 

1 Or.  E.  Rybicki 

General  Electric  Company,  Schenectady,  New  York  12010 
1 ATTN;  Hr.  A.  J.  Brothers,  Materials  and  Processes  Laboratory 

General  Electric  Company,  Schenectady,  Hew  York  12309 
1 ATTN;  Nr.  H.  F.  Bueckner,  Large  Steam  furbine  Generator 
Department 

1 Mr,  S.  Yukawa,  Metallurgy  Unit 

1 Mr.  E.  £.  Zwlcky,  Jr. 

General  Electric  Company,  Knolls  Atomic  Power  Laboratory, 

P,  0.  Box  1072,  Schenectady,  New  York  12301 
1 ATTN:  Mr.  F.  J.  Mehringer 

1 Mr.  L.  F.  Coffin,  Room  K41-K1,  General  Electric  Corp., 

P.  0,  Sox  8,  Schenectady,  New  York  12301 

United  States  Steel  Corporation,  Monroeville,  Pennsylvania 
Pennsylvania  15146 
1 ATTN:  Mr.  5.  R.  Novak 

1 Or.  A.  K.  Shoemaker,  Research  Laboratory,  Hail  Stop  78 

Westinghouse  Electric  Company,  Bettis  Atomic  Power  Laboratory, 

P.  0.  Box  109,  West  'llfflin,  Pennsylvania  15122 
1 ATTN:  Mr.  H.  L.  Parrish 

Westinghouse  Electric  Company,  Research  Laboratories, 

Pittsburgh,  Pennsylvania  15235 
I ATTN;  Or.  R.  J.  Bratton 

Brown  University,  Providence,  Rhode  Island  02912 
I ATTN-  Prof.  J.  R.  Rice 

1 Prof.  W.  N.  Findley.  Division  of  Engineering,  Box  0 

1 Prof.  P.  C.  Paris 

Camegle-Mellon  University,  Department  of  Mechanical 
Engineering,  Schrnley  Park,  Pittsburgh,  Pennsylvania  15213 
I ATTN:  Or.  J.  L.  Swedlow 
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1 Prof.  J.  D.  Lubihn,  Colondo  School  of  Mines,  Golden, 

Coloredo  80401 

1 Prof.  J.  Dvorek,  Chemtcil  Enjineerinp  Depertaent,  Duke 
University,  Durtiaa,  North  Cerolint  2770G 

George  Heshlngton  University,  School  of  Engineering  end 
Applied  Sciences,  Heshlngton,  D.  C.  200S2 
1 ATTN:  Dr.  H.  Llebowitz 
I Prof.  A.  M.  Freudenthel 

Lehigh  University,  Bethlehoi,  PennsylvanU  1801S 

1 ATIN:  Prof.  G.  C.  Sih 

1 Prof.  R.  Roberts 

1 Prof.  R.  P.  Hel 

1 Prof.  Richard  M.  Spriggs,  Assistant  to  the  President 

and  Professor 

1 Prof.  D.  P.  H.  HasselMii 

Massachusetts  Institute  of  Technology,  Cambridge, 

Massachusetts  02139 

1 AHN:  Prof.  B.  L.  Averbach,  Materials  Center,  13-S082 

1 Prof.  F.  A.  McCIIntock,  Room  1-304 

I Prof.  R.  M.  Pelloux 

I Prof.  T.  H.  H.  Plan,  Department  of  Aeronautics  and 

Astronautics 

1 Prof.  A S.  Argon,  Room  1-306 

1 Prof.  P..  L.  Coble 

Syracuse  University,  Dept,  of  Chamlcal  Engineering  and 
Metallurgy,  lletallurgy  - 409  Link  Hall,  Syracuse, 

New  fork  13210 
1 ATTN:  Mr.  H.  U.  Liu 

I Or.  V.  Heiss,  Metallurgical  Research  Labs.,  Building  0-6 

1 Prof.  E.  R.  Parker,  Dept,  of  Materials  Science  and  Engineering, 
University  of  California,  Berkeley,  California  94700 

1 Prof.  H.  Goldsmith,  Dept,  of  Mech.  Eng.,  University  of 
California,  Berkeley,  California  94720 

I Prof.  A.  J.  McEvlly,  Metallurgy  Oept.  U-136,  University  of 
Connecticut,  Storrs,  Connecticut  06268 

1 Prof.  0.  Orucker,  Dean  of  School  of  Engineering,  University  of 
Illinois,  Champaign,  Illinois  61320 

University  of  Illinois,  Urbans,  Illinois  61801 
I ATTN:  Prof.  H.  T.  Corten,  Oept.  of  Theoretical  and  Applied 
Mechanics,  212  Talbot  Laboratory 

1 Prof.  T.  J.  Dolan,  Dept,  of  Theoretical  and  Applied  Meek 

1 Prof.  J.  Morrow,  321  Talbot  Laboratory 

1 Hr.  G.  M.  Sinclair,  Oept  of  Theoretical  and  Applied  Mech. 

1 Prof.  R.  I.  Stephens,  Materials  Engineering  Olv.,  University  of 
Iowa,  Iowa  City,  Iowa  62242 

1 Prof.  0.  K.  Felbeck,  Oept.  of  Mechanical  Engineering,  University 
of  Michigan,  2046  East  Engineering,  Ann  Arbor,  Michigan  48109 


No,  of 

Copies To 

1 Or.  M.  L.  HIIIIms,  Dean  of  Engineering,  240  Benedum  Hall, 
University  of  Pittsburgh,  Pittsburgh,  Pennsylvania  15261 

1 Prof.  A.  Kobayashl,  Oept.  of  Mechanical  Engineering,  Fu-10, 
University  of  UasI  Ington,  Seattle,  Hashlngton  98195 

Albany  Metallurgy  Research  Center,  Albany,  Oregon  97321 
1 ATTN:  Mr.  R.  R.  Hells,  Research  Director 

Ford  Motor  Company,  Turbine  Research  Department,  20000  Rotunea 
Drive,  Dearborn,  Michigan  48121 
1 ATTN:  Mr.  A.  F.  McLean 

IIT  Research  Institute,  10  Hest  35  Street,  Chicago,  Illinois  60616 
1 ATTN:  Mr.  S.  Bortz,  Ceramics  Research 

United  Aircraft  Research  Laboratories,  Eust  Hartford, 

Connecticut  06108 
1 ATTN:  Dr.  Michael  DeCrescentI 

SUte  University  of  New  Fork-,  College  of  Ceramics  at  Alfred 
University,  Alfred,  New  York  14802 
1 ATTN:  fir.  Satyavan  Shukla,  Assistant  Librarian 

SKF  Industries,  Inc.,  Engineering  and  Research  Center, 

1100  1st  Avenue,  King  of  Prussia,  Pennsylvania  19406 
1 ATTN:  Harlsh  Dalai 

Deposits  and  Composites,  Inc.,  1821  Michael  Drive, 

Reston,  Virginia  22090 
1 AHN:  Mr.  Richard  E.  Engdahl,  President 

AIResearch  Mfg.  Co.  of  Arizona,  402  South  36th  Street,  Phoenix, 
Arizona  85034 

1 ATTN:  Chief,  Materials  Engineering  Department,  Dept.  93-39M 

Monsanto  Research  Corp.,  Dayton  Laboratory,  1515  Nicholas  Road, 
Dayton,  Ohio  45407 

1 ATTN:  Mr.  Richard  J.  JanowleckI,  Research  Group  Leader 

Norton  Company,  1 New  Bond  Street,  Horcester,  Massachusetts  01606 
1 ATTN:  Or.  H.  R.  Baungartner 

Federal  Mogul  Corporation,  Bearing  Group  Research,  3980  Research 
Park  Drive,  Ann  Arbor,  Michigan  48104 
1 AHN;  Hr,  P.  E.  Cowley 

Energy  Research  and  Development  Adalnlstratlon,  Division  of 
Transportation,  Highway  Vehicles  Systems  Brareh,  20  Massachusetts 
Avenue,  NH,  Hashlngton,  0.  C. 

1 ATTN:  Mr.  G.  Thur 
1 Mr.  R.  Schulz 

State  University  of  New  York  at  Stony  Brook,  Stony  Brook, 

New  York  11790 

1 ATTN:  Prof.  Fu-Pen  Chaing,  Dept,  of  Mechanics 

Director,  Army  Materials  and  Mechanics  Research  Center, 

Hatertown,  Massachusetts  02172 

2 AHN:  ORKNR-PL 

1 ORKMR-AG 

1 .Author 
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